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Abstract 

A linear shift- invariant preprocessing technique is described which re- 
quires no specific knowledge of the image parameters and which is sufficiently 
general to allow the effective radius of the composite imaging system to be 
minimized while constraining other system parameters to remain within speci- 



Introduction . The finite size of the point-spread function of a practical 


imaging system causes the output to be a two-dimensional spatial smearing or 
blurring of the original image. Methods of preprocessing the image to reduce 
this blurring have been widely studied but no fully satisfactory processor has 
as yet been developed. The preprocessing method proposed here, although new 
in formulation, can be thought of as a generalization and extension of previous 
work by Smith [4] and Stuller [ 5 ]. This method accomodates the problem of non- 
circular ly symmetric imaging system point-spread functions, provides for con- 
trolled extent of the preprocessing filter to minimize distortion due to tran- 
sients resulting from truncation errors and edge effects, can be used with 
various kinds of system noise, and can be readily extended to provide constraint 
of other system parameters. The analysis presented here relates to a line- 
scanner system although it is applicable in principle to many other system 
configurations . 

Analysis of Preprocessing Technique . Referring to Fig. 1 and noting that a bar 
over a variable denotes a two-dimensional spatial vector, a basis for choosing 
the preprocessing filter point-spread function, h r (v), in order to increase the 
resolution of the composite imaging system is to minimize the radius of gyration, 
R r , of the composite imaging system point-spread function, g(v) ; given by 
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I |v| g(v) dv 

R r 2 " ( 1 ) 

r r P - - 

J g (v) dv 
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where the composite system point-spread function is expressed 

g(v) = J h r (z) h b (v-z) dz . ( 2 ) 
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The magnitude of R f may be minimized by setting the denominator equal to a con- 

_ o 

stant and minimizing the numerator with respect to h^C'). The function |v| in 
(l) is a penalty weighting function which causes h^Cv) to decrease as v increases. 
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Additional flexibility in controlling the shape of g(v) is obtained by using 
a more general penalty weighting function w(v) which increasesat an arbitrarily 
controlled rate outside a specified range of v. With this penalty function the 
numerator and denominator of (l) can be written as the functionals 

I o “ J s 2 (v) * 6) 

-CD 

X 1 “ J g2 ^ dv = •- 00 

Two additional constraints are considered essential for most imaging sys- 
tems. The first limits the size, and to some extent the shape, of the prepro- 
cessing filter point-spread function, h^Cv), so that truncation errors are mini- 
mized. This constraint is 


00 

t 2 " J s(v) h r 2 (v) dv = K 2 (5) 
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where s(v) is a penalty weighting function similar in form to w(v). The second 
additional constraint limits the total noise power in the preprocessed image 
as follows 


.CO t .oo 

J J* 


(z-v) h (z) h (v) dz dv 
nn r 7 r ' 
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( 6 ) 


where it is assumed that the noise source is a stationary ergodic random process 
with autocorrelation function R (z-v). 

The required point-spread function is found by minimizing, with respect to 
h (•)» the augmented functional formed by adding to (3) the constraints (4), (5) 
and(6) each of which has been multiplied by an appropriate Lagrange multiplier [ 1 , 2 ]. 
The required minimization is obtained by converting to the spatial frequency 
domain, setting the gradient of the augmented functional equal to zero and solving 
the resulting equation subject to the constraints of ( 4 )- (6). The gradient of 
the augmented functional is as follows 



3 


J'TKX-v H b(5*'X H <vvy5 y ) + "X-vyV] 

-CO*CO ' + wm 

+ 'JX'WV + S *(V £ x’V f y>I} H r<?x’5y> d£ x d ?y 
+ [ax 1 |H b (f x .f y )| a + H r (f x ,£ y ) = 0 (7) 

where H^C')* W(* ), S(*), H r (*) and J(* ) are the two-dimensional Fourier trans- 
forms of h (•), w(*), s(*)» h (•) and R (• ) respectively and X-, » and \ are 
Lagrange multipliers associated with constraint equations (U)-(6) respectively. 
The frequency domain representation of the constraint equations is 


*1 ’ jTlV £ x-V| 2 | H r (f x- £ y>| 2 «x df y * K 1 < 8 > 

— eo-w 

£ 2 * JTJ r* ( V W ?y )H r ( «*• ? y > H r*< £ x- £ y > d ? x d ?y d£ x d£ y = K 2 ( 9 ) 


00 00 00 00 


— 00 — 00 — 00—00 


to 00 


I. - f [ T (f ,f ) Jh (f ,f )r df df = k . 
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y' • r * x' y ' 1 x y 3 
In order to proceed further an appropriate form for the penalty functions 
w(v) and s (v ) must be chosen. Several criteria must be satisfied by these 
functions: they must increase rapidly away from the origin; they must be Fourier 
trans f ormable ; and w(v) must be convex to assure the existence of a global 
minimum of (3). A convenient representation centered at the origin is 


2k 2k 

W(v) = W x (x)W y (y) ={(^-) “*+ “ y + c„ y } 


( 11 ) 


The constants 0 < c ,c <1 are introduced to assure that the overall weighting 

wx wy 

function will be everywhere greater than zero thus avoiding singular behavior 

of the solution. The constants k and k are positive integers. A similar 

wx wy 

expression may be written for s(v) = s (x)s (y). The shape of w(v) along one 

x y 

dimension is shown in Fig. 2 . For the following analysis it is assumed that 

k , k , k . and k are all unity. This assumption does not limit the general 
wx wy sx sy 

applicability of (7)-(l0) but is made to illustrate the method of solution of 


these equations. 
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With the above assumption (7) can be put into the following form [3] 
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( 12 ) 


Equation (l2) is a two-dimensional fourth-order partial differential equation 
which, when solved simultaneously with (8), (9), and (lO), specifies the spatial 
frequency spectrum of the required preprocessing filter. 

Image Collecting Systems with Separable Point-Spread Functions . For many 
physically realizable image collecting systems, certain assumptions may be made 
which greatly simplify the form of (12). A large class of physically realizable 
image collecting system point-spread functions may be modelled as separable 
point-spread functions, from which it follows that 

H, (f ,f ) = H, (f ) H. (f ) (13) 

b' x y bx' x 7 by y J 7 

H (f ,f ) = H (f ) H (f ) (14) 

r'xy rxxry'y 7 7 

I (f ,f ) - J (f ) J (f ) . (15) 

inn x y J-nnx x inny y ' ' 

With these assumptions, it can be shown that the solution for the two- 

dimensional point-spread function of h^*) defined by (12) and constraint equa- 

tions (8)- (lO) reduces to two sets of equations [3]: one set specifying the 

x-axis component of h^O); the other the y-axis component. Since these equa- 
tion sets are identical except for appropriate changes in x- and y-axis component 


variables, only the equation set for the x-axis components will be given. 
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This equation set is composed of the second-order ordinary differential equation 

2 H* (f ) H/ (f )x 2 

H' (f ) + b * * 0 bX 0 * * 5 H' (£ ) 
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and the constraint equations 


J ! H bx( £ x’| 2 | H rJ £ I[ )| 2 d£ x * K l, 
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- - \ -- 1 r*H^ (f ) H (f ) + H". (f ) H (f ) I df 
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+ c r |H (f )r df = K 
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f I (f ) |H (f )| 2 df = K 
J inn' x' 1 r x' 1 x 3x 

“00 

( 19 ) 

[ (f ) = H (f ) + J H (f ) . 

rx x rrx x J rix x 

(20) 


where 


Thus the spatial frequency transform of the x-axis component of h r (- ) is 
obtained by the simultaneous solution of (l6) and the constraints of (l7)-(l9). 
This system of equations can be solved by standard numerical procedures and 
the functions h^C*) and h (* ) determined. The preprocessing operation thus 
becomes two one- dimensional convolutions carried out sequentially rather than 
a single two-dimensional convolution as would generally be the case. 

Conclusions . The proposed method for specifying an image preprocessing point- 
spread function is sufficiently general to allow a variety of penalty weighting 
functions to be used for control of the effective radius of the point-spread 
function while constraining other system parameters to remain within specified 


limits. 
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Although (12) in conjunction with (8)-(l0) specify the general form of the 
spatial frequency transform of the optimum preprocessing filter, it was shown 
that in the case of a separable point- spread function the solution can be con- 
siderably simplified. This same simplification also results for the case of 
a radially symmetric point-spread function, however, the resulting processor 
must perform a two-dimensional convolution in this case rather than two sequential 
one- dimensional convolutions. 

The simplified solution based upon separable point-spread functions can 
be extended to a large class of practical systems by approximation techniques. 

For example, by properly selecting the parameters of a two-dimensional Gaussian 
function, an elliptical or circular point-spread function can be approximated 
and preprocessing carried out as sequential one-dimensional convolutions; thus 
substantially reducing the processing time. 
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Figure 2. One-Dimensional Penalty Weighting Function. 
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